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We investigate the QCD effects in the production of neutral Higgs bosons via bottom quark 
fusion in both the standard model and the minimal supersymmetric standard model at the CERN 
Large Hadron Collider. We include the next-to-leading order (NLO) QCD corrections (including 
supersymmetric QCD) and the threshold resummation effects. We use the soft-collinear effective 
theory to resum the large logarithms near threshold from soft gluon emission. Our results show 
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that the resummation effects can enhance the total cross sections by about 5% compared with the 



(N ■ NLO results. 
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I. INTRODUCTION 



The understanding of Electroweak Symmetry Breaking (EWSB) plays a key role in cur- 
rent research of high energy physics. In the Standard Model (SM), a complex scalar doublet 
is responsible for the generation of gauge bosons and fermions masses by the Higgs mech- 
anism. One neutral Higgs boson (h) survives after EWSB, which is the last elementary 
particle yet to be found in the SM. Direct searches at LEP2 set a lower bound on the SM 
Higgs boson mass > 114.4 GeV (at 95%CL) [1|, while electroweak precision measure- 
ments prefer a light Higgs boson of m h < 180 GeV [2|. 

In the most popular extensions of the SM, e.g., the Minimal Supersymmetric Standard 
Model (MSSM), two Higgs doublets are required in order to preserve supersymmetry (SUSY) 
and anomaly cancellation. In the MSSM, the Higgs sector consists of five physical Higgs 
bosons: the neutral CP-even ones h and H, the neutral CP-odd one A, and the charged 
ones H ^ . The lightest one h behaves like the SM one in the decoupling limit 3> M^o). 
Its mass is constrained by a theoretical upper bound of Mh < (130 — 140) GeV when taking 
into account the radiative corrections jjj. At lowest order, two parameters are required to 
describe the MSSM Higgs sector, which are generally chosen to be and tan/3 = v 2 /vi, 
the ratio of the two vacuum expectation values. 

For large value of tan /3, the bottom-Higgs Yukawa coupling can be considerably enhanced, 
thus Higgs production associated with bottom quarks may be quite important in MSSM. 
There are two approaches for calculating cross sections involving bottom quarks: the four 
flavor number scheme (4FNS) and five flavor number scheme (5FNS) [4]- In the 4FNS, there 
are no initial state bottom quarks. Due to non zero of bottom quark mass, large logarithms 
may appear from gluon splitting, hence transverse momentum (jpx) and pseudo rapidity cuts 
on final state bottom quarks are needed to eliminate these large logarithms. In the 5FNS, 
initial state bottom quarks are treated as massless, and by introducing a perturbatively 
defined bottom quark parton distribution function (PDF), large logarithms are resummed 
through the DGLAP evolution equation. Except for qq — > bbh, the other relevant production 
mechanisms depend on the final state being observed Q]. For inclusive Higgs production 
without bottom tagged in final state, the lowest order process is gg — > (bb)h (bb — > h) in 
4FNS (5FNS) 6|. However, if at least one high-p^ b quark is required to be observed, the 
leading partonic process is gg — > b(b)h {gb — > bh) in 4FNS (5FNS) [7J, and if two high-p^ 
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b quarks are required, the leading subprocess is gg — > bbh j§| and can only be calculated in 
4FNS. There are extensive comparison between 4FNS and 5FNS and good agreement has 



a. 



been found between the two schemes within theoretical uncertainties 

In recent years much effort has been made to the precise prediction of the inclusive 
cross section for Higgs boson production via bottom quark fusion, with neither bottom 
quark detected bb — > h. The next-to-leading order (NLO) QCD corrections Q, Q, Q and 
the next-to-next-to-leading order (NNLO) QCD corrections |l2|] to this process have been 



calculated. Also the SUSY QCD and SUSY electroweak corrections to this process have 

n 

been studied |13l |. 

When the hard scattering process involves two very different scale, the fixed order per- 
turbation expansion contains large logarithms of scale ratio. These terms might s poi l the 
reliability of the perturbation expansion and need to be resummed to all orders [l4j, [15]. 



In general, the large logarithms mentioned above can appear when the Higgs boson in the 
final state has small transverse momentum or is produced near threshold. The transverse 
momentum resummation were calculated in Refs. , and the threshold resummation 



effects are considered 
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19j with the conventional method Jj], [l5], where partial next 



next-next-to-leading-order (NNNLO) results are obtained by expanding the resummed cross 
sections. In this paper, we will further study the complete next-to-leading-logarthmic (NLL) 
threshold resummation effects on the production cross sections using soft-collinear effective 
theory (SCET) 20] in both the SM and MSSM at the CERN Large Hadron Collider (LHC). 

The paper is organized as follows: In Sec. II and III we present the analytic results at fixed 
order. In Sec. IV we use SCET to derive the resummed formula for the cross sections. In 
Sec. V the numerical results are presented and discussed. Sec. VI contains a brief summary 
and conclusions. 



II. THE LEADING ORDER RESULTS 

We consider the inclusive process A{p a ) + B{p h ) — > H^q) + X, where A and B are the 
incoming hadrons with momenta p a and pt,, = {h, H, A} are neutral CP-even or CP-odd 
Higgs bosons, with momentum q, and X is arbitrary hadronic state. 

At hadron colliders the total cross sections can be factorized into the convolution of the 
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partonic cross sections with appropriate PDFs: 

°( M i, S ) = ^2 dXa dXb fa/A(Xa,^f)fb/B(Xb,fif)^ab(z,Mf\a s (^),l4,fM 2 f ), (1) 

0,6 Jr l x " 

where a a b is the cross section for the partonic subprocess a(p a ) + b{j>b) — ► + Pa and 
are the momentum of the incoming partons a, h. The momentum fractions x a and are 
defined by x a ^) = p a (b)/Pa(b), Mi is the mass of the final state Higgs boson, r = Mf/s, and 
the scaling variable z = M?/s, where s = (p a +pb) 2 - f P /H(x,fif) is the parton distribution 
function which describes the probability of finding a parton p with momentum fraction x 
inside the hadron H at factorization scale fif. The sum is over all possible initial partons. 



Hi 




FIG. 1: Leading-order Feynman diagrams for bb — ► Hi 

The leading-order (LO) Feynman diagrams is shown in Fig. [H and its LO amplitude in 
n = 4 — 2e dimension is 

iM B = V £ rYbKiv(pb)r ii} u(p a ), (2) 

where K& = —igrrib / (2mw) is the SM coupling of Higgs boson to bottom quark, and // r is a 
mass parameter introduced to keep the coupling constant g dimensionless. Furthermore, we 
have T^ h ' HS) = 1 for scalar (h, H) production and = for pseudoscalar (A) production. 
The explicit expressions for 

sin a cos a 

n h = -, k h = -, K A = -ltWl(3. 

cos p cos p 

Here a is the mixing angle between the weak and the mass eigenstates of the neutral CP-even 
Higgs boson sector. 

The LO partonic cross sections are given by 



GB 2s 



r § J £l>k| a ^fiW = *o*(i-s), (3) 



where indicates the summation over final states and the average over initial states, 
J dPS^ represents the phase space integration and <j = ii 2e 'K\Y b Ki\ 2 / {Qs). 



III. NEXT-TO-LEADING ORDER CALCULATIONS 



The NLO QCD and SUSY QCD correctons to this process have been studied in 
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an 



131 . 127j , but we recalculate it here to check the relevant results in the previous literatures 
and make our paper self-contained. At NLO, the QCD and SUSY QCD corrections consist 
of the following contributions: the exchange of virtual gluon or gluino and the corresponding 
renormalization counterterms, the real gluon emission subprocesses, the gluon initiated sub- 
processes, and the contributions of Altarelli-Parisi(A-P) splitting functions. In the following, 
we will calculate these contributions separately. We use dimensional regularization (DREG) 
in n = 4 — 2e dimensions to regulate all divergences, and adopt MS renormalization and 
factorization scheme to remove the ultraviolet (UV) and infrared (IR) (including soft and 
collinear) divergences of QCD corrections, while SUSY QCD corrections are renormalized 



in on- shell scheme 



2l|. 



A. Virtual corrections 



The amplitude corresponding to virtual gluons exchange is given by 



M 



QCD 
V 



( 47T/i 



T(l-e) 



2 3 

+ 



(4) 



'4tt \Mf J r(l-2e) V em euv em ' 3 
where Cp = 4/3 is the SU(3) color factor. The above amplitude contains both UV and IR 
divergence. The renormalized QCD amplitude can be written as 



X QCD 



£>>■ 



with 



22 ] 



<K QCD n « s (4tt)T(1 - e) 3 
— -Of 



m b 4tt T(l — 2e) euv 

The amplitude corresponding to virtual gluino and squark exchange is given by 

M S ™ Y = M B F®. 

Here F® (i = h, H, A) are the SUSY form factors: 



(5) 



(6) 



(7) 



p(h,H) 



Of. 



C F m~ g [sin 29- b (g$ H) C (1, 1) - G%' H) C (2, 2] 



Y b K htH 4tt 

+ cos2^ (G& H) C (1,2) + G% H) Co(2, 1) 



(8) 



= ^~C F m- a (g^C (1,2) - G%>C (2, 1)) 
where Co(i,j) are the usual Passarino-Veltman three-point function 23] 



(9) 



(10) 



iGj^ are the couplings between Higgs boson and sbottom mass eigenstates, which are given 
in the appendix, 2 are the sbottom masses, m g is the gluino mass, and R b is a 2 x 2 
matrix defined to rotate the sbottom current eigenstates into the mass eigenstates: 









|. 


\b2) 









cos #5 sin 0£ 
- sin 9i cos 6^ 



(11) 



with < &i < 7r by convention. Correspondingly, the mass eigenvalues and (with 



m 



bl < m b 2 ) are § iven b y 



m~ 

bi 



mf L X b m b 



with 



= R b Ml{R b )\ Ml = 
mf I b b \X b m b m? 

02 / \ Or 



m| = M| + m 2 b + m| cos 2/3C bL , 



(12) 



m? fl = M 2 D + m 2 b - m 2 z cos 2/5C 6i? , 
X 6 = A 6 -//tan/3. (13) 

where CVl = — 1/2 + sin 2 6*^/3, CVr = sin 2 6^/3, and M| is the sbottom mass matrix. Mqjj 
are soft SUSY-breaking parameters, A b is the trilinear Higgs-sbottom coupling, and /i is the 
Higgsino mass parameter. 

The renormalized SUSY QCD amplitude is 



and the 



-My USY = A4 USY + ( + l(SZ bL + 5Z bR ) ) M B , (14) 
renormalization constants in the on-shell scheme are fixed to be 



6mf VS Y 
m b 



I i=i 



771 ~ 

B 1 9 -sm26i(-iyBo 

m b 



I 2 2 2 > 

{m hl m- gl m lit 



2 2 2 



^L = ^G F ^(4) 2 ^iK,^, 

i=l 
2 

5Z 6Ji = ^^^(^^(m^m^ml), 



where Bq } i are the two-point integrals [23]. 

After adding the counterterms, the UV divergences in A4y CD + A4y JSY are canceled, but 
the IR divergent terms still persist. The partonic subprocess cross section is 



oy =a 5(l — z) < 1 + 



r — 

F 2n \M? J r(l-2e) 



2 +2- 



5m? CD 



(15) 



+ A 



SUSY 



"ir e uv em 6 j rai, 

where A SUSY is the contribution from SUSY QCD corrections only, which is free of diver- 
gences 



A 



SUSY 



Jm SUSY 
) g 



+ 5Z bL + 5Z bR + 2F^. 



(16) 



B. Real gluon emission and gluon initiated subprocesses 



The partonic cross section of real gluon bremsstrahlung are 

a s ( 47r/i^\ e T(l — e) 



&r =a C F: 



2tt V M? / 1+1 - 2e) 



—5{l-z) 

: IR 



+4(1 + z 1 



ln(l - z) 

1 - z 



The "plus" function in Eq. (fT7|) is defined as 



dz f+(z)g(z) 



em (1 - 2) 
2 I ) In : + 2(1 - .:; 



dzf(z) (g(z)-g(l)) 



(17) 



(18) 



where g(z) is any well-behaved function in the region < z < 1. 

Combining the contributions of the LO result, the virtual corrections and the real gluon 
bremsstrahlung, we obtain the bare NLO partonic cross section: 



bare 
bb 



0+ + CTR 



a <j (1 + A SUSY )5(1 -z) + C F ^ (4tt 



rn 



2^ 



r(l - 2e) 



e m V (1 - z) + 2 J \3 



2(1 -z) 



-2 In 



1 + z 2 
M?(l-z) 



/ln(l — z) 

z 



_ 2 |l±^llnz 



(19) 



Now the soft divergences coming from virtual gluons and bremsstrahlung contributions have 
canceled exactly according to the Bloch-Nordsieck theorem 24| . The remaining divergences 
are collinear. 
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In addition to the real gluon bremsstrahlung subprocess, there are also contributions from 
the gluon initiated processes, which can be written as 

1 T(l-e) 



"bare 



e m r(l-2e) 



+ ln 



m? a - z) 



i(l-z)(7z-3) 



(20) 



The bare partonic cross sections in Eqs. f|T9|) and (120]) . which contain the collinear singu- 
larities generated by the radiation of gluons and massless quarks, have a universal structure, 
and can be factorized into the following form to all orders of perturbation theory: 

a^ re (z, I/em) = ^ T ca(z, fif, l/em) ® T db (z, ///, l/em) <g> a cd (z, ///) , (21) 

where /x/ is the factorization scale and <8> is the convolution symbol defined as 

- 1 dy 



f(z)®g(z) 



-f(y)g x 
y \y 



(22) 



The universal splitting functions T cd (z,{if, l/em) represent the probability of finding a par- 
ton c with fraction z of the longitudinal momentum inside the parent parton d at the scale 
fif. They contain the collinear divergences, and can be absorbed into the redefinition of the 
PDF according to mass factorization 25| . Adopting the MS mass factorization scheme, we 
have to 0(a s ) 



rcd(z,fif, 1/ern) = <W5(1 



1 a s r(l - e) / 4tt^ 



e m 2vrr(l-2e) 



(23) 



where Pj: d \z) are the leading order Altarelli-Parisi splitting functions [26|: 



[(1 - zf + z-] . 



(24) 



After absorbing the splitting functions T cd (z,fif, l/em) into the redefinition of the PDFs 
through the mass factorization, we derive the hard scattering cross sections a a b(z, fif), which 
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are free of collinear divergences, and depend on the scale [if. 



%0, /v, Hf) = £o < (1 + A SUSY )5(1 - + 



2tt 



2P£>(*)ln 



M 2 



// / 



+ff F *(l - z) (31n ^ + ^ 2 - 2 ) + AC F (1 + z>) (^^) 



1 

-2C F -±— In z + 2C F (1 - z) 



l-z 



(25) 



a bg {Z,fJi r ,flf) = (T 



2^ 



'/'/ 



(26) 



Finally, we combine these finite a ab (z, /x r , ///) with the appropriate partonic distribution 
function to arrive at the NLO cross sections: 



a 



NLO 



fb/A(x a ,fif)f- b/B (x b ,fif) + (X a «-> X 6 ) 



a bi (z,fi r ,fi f ) 



+ 



+ 



fb/A{x a , Vf)fg/B(x b , flf) + (x a <-> X 6 ) 
fb/A(Xa,Hf)fg/B(Xb,Hf) + (x tt «-> X 6 ) 



<^fec/(^; AV, A 1 /) 



(27) 



This result have been obtained before and our result agrees with those in Refs 



13 



27] . 



.ay 



12, 



IV. THRESHOLD RESUMMATION 



The NLO results contain terms like [ln(l — z)/(l — z)] + and 1/(1 — z) + , which are large 
near the "partonic threshold region" z — *■ 1. Physically, these singular terms represent 
a class of large logarithms of scale ratios, which come from the incomplete cancellation 
between real gluon emission and virtual gluon corrections. These logarithms can be system- 
aticall y re summed to all orders by solving the certain evolution equations in Mellin moment 



space 



14 



151 ] . One drawback in the traditional resummation formalism is that the sepa- 



ration of the contributions from the different scales is not obvious, and some ingredient in 
the resummed exponent is not easily identified with a field-theoretical object. In SCET, the 
resummation procedure has a more transparent meaning. Once the factorization properties 
is established, a soft scale of interest is separated from the underlying hard scale in the 
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framework of effective theory. By evolving from the hard scale to the soft scale through the 
renormalization group (RG) equation, the large logarithms can be resummed to all orders. 
In this approach, all the ingredients needed have a clear effective field theory interpretation. 



In fact, resummation in SCET have been carried out in deep-inelastic scattering 
Drell-Yan production [30, [3JJ, Higgs production 
and heavy colored particle production [341 ]. 



28 



32], thrust rate in e + e annihilation 



The starting point in effective theory approach to threshold resummation is the factor- 



ization formula for hadronic cross section 
<t(M?,s) = a 



30 



35 



36|, 



Ef 1 dx a f 1 d,Xb £ , \ e I \ 

/ / Ja/A{X a ,fi f )f b/ B{X b ,fif) 

ab Jt X a Jr/xa X b 



x\C v (M?,l*f)\ 3 S(VS(l-z),Mi 



(28) 



where o"o = 7r|Y;,/T;j| 2 /(6s) is the LO total cross section, Cy is the Wilson coefficient of 
operator in SCET and S is soft function. The convolution formula in Eq. (1281) can be 
further transformed into product formalism with Mellin transformation 

a N {M^s) = [ rfrr Af -V(M 2 ,s) 

(29) 



°o Yl $a N)fb/B ( x b, Hf)\Cv (M- ,ii f )\ 2 S N {M?,n f ), 

a,b 



As mentioned above, Eq. ( 1281) and (1291) contain large logarithms near partonic threshold, 
which need to be resummed to all orders. In the following, we will derive the evolution 
equations for the hard matching coefficient Cy and soft function S, respectively, in order to 
resum these large logarithms. Before proceeding, it should be pointed out that in principle, 
the threshold resummation for Higgs production through bb fusion can be obtained from 
similar results for Drell-Yan production 30|, [3J| , by replacing the hard matching coefficient 
Cy with Eq. ( 1321) . This is due to the fact that the IR divergences in our case do not depend 
on the explicit structure of the vertex, and SUSY QCD corrections do not give rise to new 
IR divergences. Nevertheless we present our full results below. 

In the full theory, the neutral Higgs boson Hi production via bottom quark fusion is 
described by the Yukawa coupling 



J{x) = Y b K4{x)Y { - i) tj){x) 



(30) 



where ip denotes the quark field coupled with the Higgs boson. In SCET this coupling can 
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be written as an effective operator 



J = C v (M 2 ,fi)Y b K^ n W n Y n T^W^ n , 



(31) 



where £ n (n) is the hard-collinear (anti-hard-collinear) bottom quark field and W n (n) (^n(n)) 
denotes the usual collinear (soft) Wilson lines which are required to ensure collinear (soft) 
gauge invariance. n and n are two light-cone vectors satisfying n 2 = n 2 = and n ■ 
n = 2. Cy{M 2 , fi) is the hard matching coefficient which comes from integrating out hard 
modes in matching from QCD to SCET. At the tree level we have Cy = 1. We can 
determine the 0(a s ) matching coefficient by evaluating the difference of on-shell matrix 
elements of operators in full theory and SCET. In the dimension regularization, the facts 
that IR structure of the full theory and SCET is identical and the on-shell integrals are 
scaleless and vanish in SCET imply that the UV divergences of SCET is just the negative of 
the IR divergences of the full theory. Furthermore, the hard matching coefficient is simply 
the finite part of the full theory virtual amplitudes. From the 0(a s ) virtual corrections 
Eq. (TT5]) we can obtain the NLO matching coefficient 



C v (Mf,ii) = l + ^C F 



-In 2 



— 71 — Z 



I _ A SUSY 

+ 2 



(32) 



M} J ' 6' 

where, for simplicity, we redefined the 't Hooft mass as /j, 2 — > /x 2 e 7B /(47r) and 7^ is the Euler 
constant. The RG equation of the above matching coefficient can be found from the UV 
divergences in the effective theory, 



d 



with 



where 71 (/i 
dimension 



din /1 
7i(A*) = E ( 



C v (M?,ri = 'y 1 (riC v (M?,n), 



n=l 



7T 



(n) 



(33) 



(34) 



is the anomalous dimension of Cy, and A\ is the well known cusp anomalous 
371 . |38| . originates from the 1/e 2 poles in the UV divergences of the effective 



theory. Note that the anomalous dimension 71 (/x) itself contains a ln(M 2 //i 2 ) term, which 
leads to Sudakov double logarithms evolution, while the term leads to single logarithms 
evolution. We can extract A\ and A at the LO from 0(a s ) virtual corrections 



4 1} = & 



A 



(i) 



3 -c F . 

2 



(35) 
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In order to reach the NLL accuracy, we need the two loop expression of A\ 



(2) 



2 



C A \-^--\--n 



(36) 



,!8 6 
where = 3 and rif = 5. 

The soft function 5, defined as the closed Wilson loop formed from the product of the 



301 ]. describes the rea. 



soft Wilson lines in the two currents 
exchange in the soft limit. At the NLO it is given by 



39 



duon emission and virtual gluon 



S(V§(l-z),M?,» a 



5(1 -z) + 
ln(l 



2n 



C f 



+ 



6(1 
4 



In 



2 

M? 



7i~ 

~2 



In 



A 

M? 



In moment space the soft function can be written as 

0<-s(Hs 



S N (Mf^ s ) = l + 



2tt 



-C F 4 In 



2 Nfi, 



6 



-7T 



(37) 



(38) 



where N = Ne lE . It is manifest in Eq. (|38l) that the contribution from the logarithms 
In TV can be eliminated by choosing the scale /i s = /// ~ Mi/N. The same scale choice 
is also adopted implicitly in the traditional approach. However, such scale choice should 
be taken with caution, it will leads to a Landau pole at N ~ Mj/Aqcd when N becomes 
large. Such spurious Landau pole singularities is avoided in the SCET approach by assuming 
Mi/N ^ Aqcd [28|. In other words, the soft function in Eq. (138!) is only applicable at a 
perturbative calculable scale. In principle, nonperturbative effects might be important and 
a modeling of the soft function at nonperturbative scale is then needed 4jJ, |42j. Never- 
theless, the anomalous dimension of the soft function given below is expected to be free of 
nonperturbative corrections, hence its evolution may still provides valuable information. A 
close investigation of the soft function in SCET is beyond the scope of the present paper, 



and we refer the reader to Ref. 



41 



421 ] for the detailed discussion. 



The evolution equation of the soft function can be derived from the fact that the cross 
section in the threshold region is independent of the factorization scale. In moment space 
the soft function obeys the evolution equation 

d 



dln.fi 



S N (Mt,fj) = (-2 7 i(/i) + 2 7 2(/i))MM/,/x), 



where 72 (fi) governs the DGLAP evolution for the PDFs 26 ] 

d 



din// 



faa(x,V) = -72(0) / OT 



(39) 



(40) 
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and similarly for f~(x,[i). It can be shown 28] to all orders in perturbation theory that 



the anomalous dimension 72 (ju) is a linear function of In N 

7fM = E[(v)"K""^ +B »")]- < 41) 



n=l 



where the expanding coefficients can be obtained from two loop splitting function [43j] 

B? = 2A?, BP = - 3 -C F , 

Bf ] = 2Af\ (42) 

Note that at 0(a s ), the term Bq in 72, which leads to single logarithms evolution, coincides 
with the corresponding term in 71. It can then be seen from Eq. ( 1391) that the anomalous 
dimension of the soft function doesn't lead to single logarithms evolution. However, this is 
no longer true at the NNLO 30|. We also notice that the peculiar IniV evolution in the soft 
function anomalous dimension corresponds to a plus distribution evolution in momentum 
space. With the extra angular restrictions on the real gluon, it originates from the incomplete 
cancellation between real and virtual corrections in the threshold region. 

Combining the above results, we can write down the resummed cross section in moment 
space 

a N N LL (M?,S) = a ^/^(x a ,/i / )/ fe %(a; b , / , / )|^(M2,^)| 2 ^(M2, /i/ ) 

a.b 

= °o E fa/A^a, Vf)f b N /B (x b , H f )\C v {M?, »r)\ 2 S N (M?, ^)e h+I \ (43) 

a.b 

with 

h = -2 %,(n), h = 2 H ^72 (/i). (44) 

Now we can evaluate the integrals in Eq. (I44p using the two-loop evolution of a s in the 
MS scheme. Keeping only terms up to NLL in the exponents in Eq. (|43l) . we obtain 

h + I 2 = ha.Ngi(P a s (fi r ) \nN/n) + g 2 {^a s (^ T ) lniV/vr) 

+0(a s (a s lnN) k ), (45) 
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with 



<?i(A) 
ft.(A) 



(i) 



AoA 



[2A + (1 -2A) ln(l -2A)] 



(46) 



Ao 
4 2) 



ln(l - 2A) + 



A 3 



o 



2A + ln(l - 2A) + - ln 2 (l - 2A) 



2 [2A + ln(l-2A)]-^ln(l-2A)ln^ 



Mi 



4 (1) ;y 2 
^2Aln^ 

Ao /i/ 



.4, 



(i) 



5, 



(i) 



0o 



ln(l-2A), 



where (3q and (3\ are the first two coefficients of the QCD f3 function 44j: 



A) = ^( n CU - 2n/), Ai = ^(!7Cl - SCUn, - 3C F n/). 
The NLL cross section in moment space is then given by 

a o ^2 fa/A( x a, t*f)fb/ B ( x b, Hf) \C v (Mf, fj, r ) \ 2 S N (M?, ///) exp(#i \nN + g 2 ) 



NLL 



(47) 



(48) 



(49) 



u.b 



Note that the NLL resummed cross section for SM Higgs boson production can be obtained 
from Eq. (|4"9"|) by setting Ki — 1 and A SUSY = 0. To obtain the physical cross section, we 
perform the inverse Mellin transformation back to the x-space 

1 



a NLL M 



2ni 



dNr~ N af LL 



'N 



(50) 



c 



Here the integra 
duced in Ref. 



contour is chosen as the minimal prescription 45j and the tricks intro- 
46] is used to evaluate the iV-integral numerically. Note that in the minimal 
prescription, the Landau pole singularities in the numerical inverse Mellin transformation 
is avoided by choosing the contour that doesn't include this pole. Another approach devel- 
oped recently is to solve the RG equation in momentum space directly 30J, which we do 
not consider here. Finally, The resummed cross section at NLL accuracy is defined to be 
the NLL cross section plus the remaining terms in the NLO result which are not resummed, 
i.e., 



RES 



NLL , NLO 
o + cr 



a 



NLL 



a s =0 



da 



NLL 



da. 



a s =0 



V. NUMERICAL RESULTS AND DISCUSSION 



In this section, we present the numerical results for inclusive production cross section 
of neutral Higgs bosons at the LHC. In our numerical calculations the following SM input 

14 



parameters were chosen 471 ] : 

M t = 172.4 GeV, G F = 1.16637 x 10~ 5 GeV" 2 , a s (M z ) = 0.1176, 
M w = 80.398 GeV, M z = 91.1876 GeV. 

The running QCD coupling a s was evaluated at the two- loop level 



(51) 



PDFs 



48j, and the CTEQ6.6M 



491 ] were used to calculate the various cross sections. Moreover, in order to im- 



prove the perturbative calculations, 1-loop and 2-loop running masses mb(fi r ) are taken as 
following [22, lad]: 



m 6 (/i, 



r J 1— loop 



m 



pole 



co /bo 



a s (ml ole ) 



(52) 



for LO cross sections and 



m b (/j, r ) 2-ioop = ml 



pole 



co /bo 



a s (m p b ole ) 



1 + -T~(ci - &i) (oc a (pr) ~ a a (m; 



poleN 
6 / 



(53) 



(54) 



for NLO cross sections, respectively, where 

h 1 f 11 A r 2 Ar \ t 1 / 5jJVc-lgiVj 

6 = 4^ iT^c - 3^/) > 6 1 = 2^ lljy e -2JV; 

c = ^, Cl = ^(101^-10^;. 

In addition, to resum the leading tan/? enhanced effects from SUSY QCD corrections, the 
LO cross sections <To is replaced by the "Improved Born Approximation" (IBA) 13j, |5CJ, [51] 
as following: 



(Jo. 



a 



IBA 



:i + a,) 2 



1 - A,, 



1 



a 



H 

IBA 



<7, 



ii 



with 



A 

a IBA 



A,, 



:i + a,)2 



(7n 



'1 + A h ) 2 



1 + A, 



tan a tan /3 
tana 



tan j3 



A, 



tan 2 ,3 



where J (a, 6, c) is defined as 
I(a, b, c) 



2^ 



1 



C F m- g \i tan (31{m~ bi , m^ 2 , m 5 ) 



j 2 

> />./• awl 2T( a ^ ln l2" + ^c 2 ln— + cVln— ). 



ft 2 



C 



(55) 



(56) 



(57) 



To avoid double counting, it is necessary to subtract the corresponding SUSY QCD correc- 
tions from the renormalization constant drrib in the following numerical calculations. 
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All the MSSM parameters are generated with FeynHiggs 52]. For simplicity, we only 



present numerical results for the m™ ax scenario, which is suitable for the MSSM Higgs 



boson search at hadron colliders 



531 ]. and the resummation effects on total cross sections for 



other scenarios are almost the same. In the m™ ax scenario [541 ]. The parameters are: 

M SUSY = 1 TeV, fj, = 200 GeV, M 2 = 200 GeV, 

X 4 = 2M SUS Y, A b = A u m3 = 0.8M SUS Y, (58) 

where X t = A t — fi cot (3 and M 2 is the wino mass term. 

Moreover, to show the resummation effects in SM (MSSM) Higgs boson production, we 
define 

SK = -,-NLO » (59) 



(J 1 

where cr NLO includes NLO QCD and SUSY QCD corrections. 

In Fig. [2] we show the K factor, which is defined as the ratios between the cross sections 
at the higher orders and the cross sections at the LO, for SM Higgs production cross sec- 
tions with the NLL resummation effects and the NNNLO collinear and soft gluon effects in 
Ref. [l9j, respectively, assuming u~ = jif = M h . We can see that our result is about 5%, 



while their result shown in Ref. [19J is about 35%. 

Fig. [3] shows the total cross sections for SM Higgs production as functions of the mass 
of Higgs boson including higher order QCD effects, assuming \i r = fif — M h . It can be 
seen from the figure that the total cross sections become small as the Higgs boson mass 
increases, which is due to the decreasing of the bottom quark density. The figure shows that 
the threshold resummation effects reduce the LO results significantly, and enhance the NLO 
results by a few percent generally. 

Fig. H] shows the renormalization scale dependence of SM Higgs production, assuming 
lif — Mh/4. We find that the renormalization scale dependence is reduced by the resum- 
mation effects, and the resummed cross section is very close to the NLO cross section in 
the vicinity of \x r ~ Mh- If restricting the renormalization scale to a factor of five above or 
below M h , the scale dependence is reduced from about 50% at LO to 38% at NLO, and to 
about 30% at NLL. 

Fig. [5] gives the factorization scale dependence of SM Higgs production, assuming \i T = 
M h . We see that the NLO corrections reduces the scale dependence significantly, and the 
resummation effects can not improve the scale dependence. This is due to the fact that 
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12] can further reduce the factorization 



the dominate contribution to the reduction of factorization scale dependence comes from 
the splitting processes of the initial states at the NLO, which are not included in the NLL 
resummation effects, while the NNLO corrections 
scale dependence. As shown in the figure, the resummation effects is small around fif ~ 
which implies that the convergence of the resummed logarithmic terms is very well 
around such scale. 

In Fig. [6] the resummation effects 5K are presented as a function of the SM Higgs boson 
mass M h , assuming /i r = M h for jif = M h /4 and fif = M h , respectively. The results show 
that the resummation effects are quite small for fif = Mh/4, about -1%, but for fif = Mh, 
SK can be over 5%. The resummation effects do not lead to large corrections relative to the 
NLO results, since the scaling variable r = M%/s is far away from 1 and the falling off of 
the bottom quark density is smooth. 

Figs. 7-12 show the total cross sections for h, H and A production as functions of their 
masses for the m™ ax scenario, assuming tan/3 = 5 and 30, respectively. In these figures and 
the following, the variation on Mh and Mh is obtained from varying Ma- We find that in all 
cases, the NLO SUSY QCD corrections are negligible after employing the Improved Born 
Approximation, and the threshold resummation cross sections reduce the LO cross sections 
and enhance the NLO cross sections by a few percent. 

Figs.[13]and[Tl]present the resummation effects 5K as a function of Mh assuming fi r = Mh 
for tan/3 = 5 and 30, respectively. In general, 5K is about -1% when fif = Mh/4. However, 
the SK can be larger than 5% when fif = Mh. Note that the SUSY QCD corrections has 
little impact on the resummation effects 5K as shown in the figures, since the SUSY QCD 
do not give rise to new soft gluon interaction. 

VI. CONCLUSION 

In conclusion, we have calculated the QCD effects in the production of the neutral Higgs 
boson via bottom quark fusion in both the SM and the MSSM at the LHC, which include 
not only the NLO QCD and SUSY QCD corrections, but also the NLL threshold resum- 



mation effects in the framework of SCET. Similar to Drell-Yan production 30, |3l|], The 
resummation is achieved by separating the contribution from hard and soft scale into dif- 
ferent matching coefficients and then summing the large logarithms of scale ratio via RG 
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equation. This approach has the advantage that the ambiguity due to the Landau pole 
singularities is reduced. Our results of the NLO QCD and SUSY QCD corrections agree 
with the calculations reported in the previous literatures, and the resummatioin effects are 
about -1% and 5% for fif = M h /A and fif = M h , respectively. 
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APPENDIX A 



In this appendix, we collect the relevant MSSM Feynman rules as following 55l. l56| 
1. The coupling between neutral Higgs boson and sbottom H^ — b\~ b* m : 

i[R i GW(R i ) T } lm = i(GW) lm , 



with 



G (h) 



g^lMl sin(a + P)C bL + mg£j a (%f- + ^ 



M W \ a ( Ab ^l + ^ -Mf sin(a + (3)C bR + mg£f 



. g 2 mfe / -A& tan j3 — \i 



G (A) 



Zrriw I A b tan(3 + n 



2. The coupling between bottom, sbottom and gluino b — b) — g: 




■V2ig s T a kl (R^P L - R) 2 P R ) 
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k , 



I 

b <- 



-yfiigsTf^&Pn - R%P L 



9 



Here Pl(r) = (1 T 7s)/2 and T a is the SU(3) generator in fundamental representa- 
tion. 
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FIG. 2: The K factor of our results compared with those of Ref 



19j |. assuming [i r = fij = M/j 
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FIG. 3: The total cross section for pp — > h + X in SM at y/s = 14 TeV, assuming /i r = fj,f = Mh. 
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FIG. 6: The NLL resummation effects with different factorization scale for pp — > h + X in SM at 
y/s = 14 TeV, assuming /x r = M^. 




FIG. 7: The MSSM total cross section for pp —> h + X at y/s = 14 TeV in the m™ ax scenario, 
assuming fj, r = fj,f = and tan/3 = 5. 
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FIG. 8: The total cross section for pp —* H + X at y^s = 14 TeV in the m™ ax scenario, assuming 
jjL r = fj,f = Mh and tan (5 = 5. 




FIG. 9: The total cross section for pp — > A + X at y/s = 14 TeV in the m™ ax scenario, assuming 
fj, r = fj,f = Ma and tan /3 = 5. 



26 



1 1 1 1 I 1 1 1 1 I 1 1 1 1 I 1 1 1 1 I 1 1 1 1 I 1 1 1 1 I 




100 105 110 115 120 125 130 

M h (GeV) 



FIG. 10: The total cross section for pp — » h + X at = 14 TeV in the m™ ax scenario, assuming 
jjL r = fif = Mh and tan j3 = 30. 




FIG. 11: The total cross section for pp — > H + X at ^/s = 14 TeV in the m™ ax scenario, assuming 
fj, r = fj,f = Mh and tan/3 = 30. 
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FIG. 12: The total cross section for pp —* A + X at y^s = 14 TeV in the m™ ax scenario, assuming 
fj, r = fj,f = Ma and tan/? = 30. 
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FIG. 13: The comparison of NLL resummation effects at different factorization scale for — > h+X 
at yfs = 14 TeV in the m™ ax scenario, assuming fi r = Mh and tan/3 = 5. 
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FIG. 14: The comparison of NLL resummation effects at different factorization scale for pp — > h+X 
at = 14 TeV in the m™ ax scenario, assuming [x r = and tan/3 = 30. 
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